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Abstract The dynamics of cold atoms in conservative optical lattices obviously depends on the geometry 
of the lattice. But very similar lattices may lead to deeply different dynamics. In a 2D optical lattice 
with a square mesh, it is expected that the coupling between the degrees of freedom leads to chaotic 
motions. However, in somme conditions, chaos remains marginal. The aim of this paper is to understand 
the dynamical mechanisms inhibiting the appearance of chaos in such a case. As the quantum dynamics of 
a system is defined as a function of its classical dynamics - e.g. quantum chaos is defined as the quantum 
regime of a system whose classical dynamics is chaotic - we focus here on the dynamical regimes of classical 
atoms inside a well. We show that when chaos is inhibited, the motions in the two directions of space are 
frequency locked in most of the phase space, for most of the parameters of the lattice and atoms. This 
synchronization, not as strict as that of a dissipative system, is nevertheless a mechanism powerful enough 
to explain that chaos cannot appear in such conditions. 
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1 Introduction 

Optical lattices are one of the most efficient tools to ma- 
nipulate cold atoms, by tuning or adjusting parameters 
such as the mesh and height of the sites (atom confine- 
ment, atomic density), or the lattice geometry. Thus it 
is not surprising that they became a toy model in many 
fields. However, depending on the considered situation, the 
role of the interactions between atoms varies a lot. Con- 
densed matter systems and strongly correlated cold atoms 
in optical lattices offer deep similarities. The flexibility 
of the latter allowed the observation of the superfluid- 
Mott insulator quantum phase transition [T] , of the Tonks- 
Girardeau regime [2] , and more generally to the superflu- 
idity properties, including the instabilities. In these ex- 
periments, the interactions between the cold atoms play 
a major role, and require the use of a Bose-Einstein con- 
densate. In particular, instabilities are described by the 
Gross-Pitaevskii equation, through the non-linear term 
[3] . Quantum computing also required a coupling between 
qubits. Optical lattices with controlled or "on demand" 
interactions appear to be an efficient implementation of 
a Feynman's universal quantum simulator [4^, and are 
among the most promising candidates for the realization 
of a quantum computer [5] . 

On the other hand, interesting behaviors can be found 
in noninteracting systems. In many situations, including 
the one considered here, the underlying physics is that 
of a single atom, without any interaction between neigh- 
bors. The higher number of atoms simply increases the ob- 
servable signal. It is the case in statistical physics, where 



cold atoms in optical lattices, through their tunability, 
made possible the observation of the transition between 
Gaussian and power-law tail distributions, in particular 
the Tsallis distributions [6l[7] . Non-interacting cold atoms 
in optical lattices also allowed the observation of Ander- 
son localization [SllQllTO]. Such cold atoms appear also to 
be an ideal model system to study the dynamics in the 
classical and quantum limits. Both are closely related, as 
the latter is only defined as a function of the former. For 
example, quantum chaos is defined as the quantum regime 
of a system whose classical dynamics is chaotic. A good 
understanding of the classical dynamics is therefore an es- 
sential prerequisite to the study of quantum dynamics. In 
non dissipative optical lattices, both the classical and the 
quantum situations are experimentally accessible, and it is 
even possible to change quasi continuously from a regime 
to the other [11]. Moreover, the extreme flexibility of the 
optical lattices makes it possible to imagine a practically 
infinite number of configurations by varying the complex- 
ity of the lattice and the degree of coupling between the 
atoms and the lattice. Many results have been obtained 
during these last years in the field of quantum chaos [TTl 
I12j . As mentioned above, these results in the quantum 
regime follow extensive studies of the classical system [T3] . 

Most of the above works used very simple potentials, 
mainly ID. For example, chaos is obtained only with a 
periodic (or quasi-periodic) temporal forcing of the am- 
plitude or frequency of the lattice [HHS], and only the 
temporal dynamics of the individual atoms is studied. But 
recently, it appeared necessary to introduce more complex 



2 



D. Hennequin and P. Verkerk: Synchronization in non dissipative optical lattices 



potentials, in particular 2D potentials [14 . Although the 
dynamics of particles in 2D potential has been extensively 
studied in the past, it was mainly in model potentials [15[ 
[16]. Experimental optical lattices approach these models 
at best on a limited domain, at the bottom of the wells. 
But in most cases, the potential is more complex, and 
leads to a more complex and richer dynamics [TT . Under- 
standing acutely the classical dynamics of atoms in real 
potentials is important, in particular because it has signifi- 
cant consequences in the corresponding quantum systems. 

The most common approach for the study of complex 
dynamics in conservative systems is statistical, e.g. eval- 
uating the percentage of the chaotic area in the phase 
space. However, a more deterministic approach is possi- 
ble, as in dissipative systems. In a recent study, we stud- 
ied the dynamics of atoms in different 2D conservative 
optical lattices. We chose experimentally feasible lattices, 
as these studies were motivated by experiments, and we 
showed that different types of chaotic dynamics appear, 
leading to different macroscopic behaviors: we showed in 
particular that the lifetime of atoms in the lattices depend 
drastically on their dynamics [17j . 

One of the simplest experimental 2D conservative po- 
tential that we studied in [T7^ is the square lattice, result- 
ing from the interference of 2 orthogonal pairs of counter- 
propagating stationary waves. This lattice has a square 
mesh, and the two directions in space are strongly coupled. 
Therefore the dynamics is expected to be fully chaotic 
when anharmonicity is high enough, i.e. for high enough 
energy of the atoms. This fully chaotic regime is effec- 
tively observed, except when the lattice is red detuned. In 
this case, chaos disappeared almost completely, and the 
dynamics remains essentially quasiperiodic, although the 
nonlinearities remain the same. The reasons of this un- 
usual behavior was not discussed in[T3, where we focused 
on the chaotic behaviors. However, the lack of chaos where 
it is expected deserves to be studied in details, to under- 
stand what are the mechanisms inhibiting its appearance? 

In the present paper, we give some answers to this 
question by studying in detail the square optical lattice 
with red detunings. We show that at the bottom of the 
wells, the resonance frequencies in both directions are de- 
generate, but when the atom energy increases, this de- 
generacy should obviously disappear because of the an- 
harmonicity of the potential. However, we show that the 
motions in both directions remain locked to the same fre- 
quency on a large domain, following a synchronization 
mechanism close to the frequency locking process of dis- 
sipative systems. Because of the conservation of energy, 
it is not a strict frequency locking, but the quasiperiodic 
regime appears to be mainly a frequency locked periodic 
regime with small sidebands. Even when the edges of the 
wells are approached, chaos appears very marginally, in 
a regime where the frequencies remain locked. The paper 
is organized as follows: after this introduction, we discuss 
the model of the lattice, and we search analytically for pe- 
riodic solutions at the bottom of the well, where approxi- 
mations lead to a Duffing-like model. Then we discuss the 
general case, i.e. far from the bottom of the well. In this 




Figure 1. a) Layout of the laser beams, b) Spatial distribu- 
tion of the intensity in the {X, Y) space for a = 0.5. Black cor- 
responds to the minimum value (zero intensity), while white 
corresponds to the maximum. The dashed square delimits the 
elementary mesh of the lattice, and the white crosses are the 
saddle points. 

case, results are obtained mainly from numerical simula- 
tions. We show that the main frequency of the motion is 
the same in both directions, whatever the energy is, and 
we analyze the mechanisms leading to this synchroniza- 
tion. 



2 Description of the square lattice 

When cold atoms are dropped in a stationary wave, they 
undergo a force F, the potential U of which is proportional 
to the wave intensity /, and inversely proportional to the 
detuning A between the wave frequency and the atomic 
transition frequency: 

F = -VU 
C/cx 4- 

The detuning is a key parameter for the behavior of the 
atoms: they accumulate in bright sites for red detunings 
(Z\ < 0), and in dark sites for blue detunings (Z\ > 0). This 
paper is restricted to the red detuned case, i.e. atoms lo- 
cated in the bright areas. When the atoms are cooled with 
a Magneto Optical Trap (MOT), the atomic density in 
these optical lattices is small enough to neglect the colli- 
sions between atoms, and so the only source of dissipation 
is the spontaneous emission. As spontaneous emission is 
proportional to I/A'^, it is relatively easy to build conser- 
vative optical lattices. 

The lattice geometry is a decisive parameter for the 
motion of the atoms. Many geometries can be experimen- 
tally built, as a vertical stack of ring traps [TS], five-fold 
symmetric lattice |19] or quasiperiodic lattices [20]. But 
one of the simplest experimental optical lattice is the case 
of two orthogonal stationary plane waves with the same 
linear polarization. The configuration of the laser beams is 
shown on Fig. [T^. The total field is £ — cos kx + e'"^ sin ky, 
where x and y are the two space coordinates, (j) a phase, 
k = 2tt/X the wave vector and A the wavelength of the 
laser beam. The intensity can be written as: 

/ = cos^ kx + cos^ ky + 2a cos kx cos ky 
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where a = cos 4>. With the adequate normahzation, the 
potential is, for A < 0: 

U = -I 

When a = 0, the coupling between x and y disappears, 
and the problem becomes separable. In all the other cases, 
the coupling between x and y could induce complex dy- 
namics. It is easy to see that in these cases the elementary 
mesh of the potential is turned of 7r/4 as compared to the 
{x,y) axes, and thus it is natural to introduce the follow- 
ing new coordinates: 

X ^ kx + ky 
Y — ky — kx 

The intensity and the potential can now be written: 

/ = -[/= 1 + a (cos X + cosy) + cos X cos r (1) 

The minus sign is that of the red detuning {A < 0) that 
we consider here. In this conditions, the atoms accumulate 
in the bright areas: the maxima of intensity correspond to 
the minima of the potential. For blue detunings, atoms 
accumulate in the dark areas, leading to significant diff'er- 
ences in the motion [T7] . 

As an example, the following of the paper will be il- 
lustrated with numerical values and plots obtained with 
a = 0.5. The choice of this value is not restrictive, and 
these illustrations are representative of the behaviors for 
other values of a. Fig. [T|d shows the spatial distribution 
of the intensity in this case. The elementary mesh is in- 
dicated through the dashed line. Assuming a > 0, the 
potential U has its main well at the absolute minimum 
Eq = —2 (1 -|- a) at coordinates (n27r, m27r), where m and 
n are integers. It has also a relative minimum —2 (1 — a) 
in (tt -|- n27r, tt -|- r7i27r). a = is a special case because the 
absolute and relative minima have the same value. The in- 
tensity goes to zero in (tt + n27r, m27r) and (n27r, tt + m27r), 
corresponding to the maxima of the potential. Two neigh- 
boring maxima are separated by a saddle point where 
U = Et = -{1 



a ) (remind that a < 1). Et is the min- 



imum energy required for a classical atom to jump from 
one well to another one. In the quantum world, atoms 
can also tunnel through the barrier, leading to a band 
structure. However, we do not have to consider the tun- 
nelling here, because, prior to any quantum treatment, 
the classical behavior has to be deeply understood. On the 
other hand, it has been shown that, for wells deep enough, 
the tunnelling vanishes [21]. Thus atoms, the energy E of 
which is smaller than the threshold ET^ are trapped into 
one site. On the contrary, atoms with E > Et can travel 
between sites, if they move in the right direction. It is 
also important to note that these saddle points are on the 
bissectors, connecting on a straight line a main well to a 
secondary one and again to the next main well. Thus an 
atom with E > Et following this straight line does not 
meet any obstacle: bissectors are clearly escape lines. 

Inside a trap site, the energy of the atom plays the 
role of a stochastic parameter. Indeed, for low energies, 
the atoms remain located close to the bottom of the well, 



and their dynamics can be approximated by an harmonic 
motion. As the energy increases, the potential becomes 
more and more anharmonic, the nonlinearities increase, 
and the dynamics can become more and more complex 



3 The lattice around the origin: periodic 
solutions 

We study in the following the classical dynamics of cold 
atoms in an optical lattice, and more generally that of 
a classical particle in the corresponding potential. Let us 
first examine what is the motion at the bottom of a main 
well. We approximate the potential to the fourth order: 



U 



{a + l){Dx + DY) 



(2) 



with 



1 
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y2 



Dx = l- — + — 
DY 



24 



24 



Dx and Dy have the shape of the potential of a conser- 
vative Duffing oscillator. Thus the motion of the particles 
at the bottom of a well appears to follow the dynamics of 
two coupled Duffing oscillator, and so we can search for 
approximate harmonic solutions. The equations to solve 
are immediately derived from eq. 



Fx = -{a + l)[X 



X^\ XY^ 



, ( Y^\ X^Y 

= Y 



(3a) 
(3b) 
(3c) 
(3d) 



Let us first search if periodic oscillations are approximate 
solutions of these equations. We search for: 



X = Xq cos {ut + ipx ) 
Y = Yo cos {ojt + iPy) 



(4a) 
(4b) 



where X^, Yq, uj, ipx and ipY are constant. First, we look at 
the constraints that have to be satisfied for the solutions. 
Then, to fully characterize a solution, we need to check its 
stability, i.e. the behavior of the trajectories close to that 
solution. In a dissipative system, a stable periodic orbit is 
an attractor, and thus it plays a crucial role in the effective 
dynamics of the system. On the contrary, an unstable pe- 
riodic orbit is not directly observable in experiments. The 
situation is not so different in a conservative system. A pe- 
riodic orbit is stable if the behavior in its vicinity changes 
slowly with the distance, i.e. if the trajectories change con- 
tinuously from the periodic orbit to a torus. The motion in 
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the vicinity of the periodic orbit is then not very different 
from a periodic osciUation. On the contrary, if the periodic 
orbit is unstable, a small change in the initial conditions 
leads to a completely different behavior, which cannot be 
considered as a small perturbation of the initial periodic 
cycle. 

The force terms in ([3]) are developed in their Fourier 
components, dropping the higher harmonics: 

Fx = Ax cos {ujt + i)x) + Bx sin {tot + ij^x) (5a) 

Fy = Ay COS {ujt + ifjy) + By sin {uit + f/'y) (5b) 

Ax,Y and Bx,y are the corresponding Fourier compo- 
nents: 



Ax Y = — dt Fx Y COS {ujt + iPxy) 

Bx,Y ^ — i dtFx,Y si-a^iuJt + il^x.Y) 



(6a) 
(6b) 



which lead to: 
Ax ■ 



XqYq 



Bx=- 



(2 + COS (2V'y - 2i/'x)) 
sin {2ijjy — 2ipx ) 

-I n 



Ay = -yo(l + a) + -^(l + a) 



B. 



8 



(2 + cos {2iPy -2ii;x)) 
sin (2-0^ — 2^x) 



(7a) 
(7b) 

(7c) 
(7d) 



Differentiating twice the equations ^ , we obtain another 
expression for the system equations which imposes Bx = 
By = 0. There are 6 periodic solutions for Eqs ([3]): 







(8a) 


Yo = 




(8b) 


i'x = i^Y 




(8c) 


■>Px =1pY ^ 


hTT 


(8d) 


'^X = IpY ~ 


TT 

~ 2 


(8e) 


^X =^Y ^ 


TT 

^2 


(8f) 



The trivial solutions (|8ap and (|8bp are those of a par- 
ticle oscillating along one of the main directions, where 
the minimum of the motion coincides with the bottom of 
the well. As the potential is invariant when X and Y are 
exchanged, we have to study only one of these solutions. 
The antiphase solution (j8dp corresponds to a change of 
sign of Y as compared to in-phase solution ((8c|) . As the 
potential is even in X and Y , we need to study only one 
of them. For the same reasons, only one of the quadrature 
solutions (l8el) and (l8fl) has to be studied. 



Trivial solutions 

Let us first look at the Yq = trivial solution. It comes 
immediately from Eqs (j7ap and jS]) that the frequency of 
the motion is: 



1 - 



(9) 



where wq = vT+a is the frequency of the oscillations at 
the very bottom of the wells. This relation is nothing but 
the well-known dependence on the amplitude of motion, 
of the period of oscillation for a simple pendulum T = 

To (l + f ). 

Let us look in the neighboring of our trivial solution. 
For Y small, Eq. (pc|) writes: 



Y + Lu^Y- 



X^Y 



= 



(10) 



With a solution of the type of pajl , this equation becomes 

f+ (^cj2_r^(l_a)_M cos2cji^y = (11) 

This equation is that of a parametric oscillator with a 
frequency close to lu, with no damping, and excited at 
the frequency 2uj. Such an oscillator is known to diverge 
rapidly, and thus the behavior in the vicinity of the trivial 
solutions does not remain close to these solutions. Thus 
trivial solutions are unstable solutions of the system. 



In-phase and antiphase solutions 



Let us now look at the ([He]) in-phase solution. Eqs ([7a|) 
and (|7cP give two expressions for w, which are compatible 
only if the motions in the two directions have the same 
amplitude Xq = Yq, with a frequency given by: 



(12) 



This solution appears to be very similar to the trivial solu- 
tions: the particles oscillate along a straight line, namely 
one of the bissectors. Their frequency decreases as the 
motion amplitude increases. However, as the stochastic 
parameter of our potential is the energy, it is more con- 
venient to represent the evolution of the dynamics as a 
function of this parameter. The conversion between am- 
plitude and energy depends on the motion. For the present 
motion, the total energy Ei leading to a motion amplitude 
of Xq is equal to the potential energy in (Xq^Xq). Thus 
the relation between Xq and Ei is easily deduced from 
Eq. 0: 



El 



^0 P^o 



2)-(" + 4)§ (13) 



Fig. [2^ shows (solid line) the evolution of lui as a func- 
tion of the energy, for a = 0.5. In this case, the bottom 
of the main well corresponds to an energy Eq = —3, while 
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-3.0 -2.5 -2.0 -1.5 -1.0 
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Figure 2. Oscillation frequencies of the particle in the vicinity 
of the periodic orbits, in the case of the Duffing approximation. 
In (a), frequency of the periodic solution; in (b), frequency of 
the beating. The solid (resp. dashed) lines refer to the analyt- 
ical in-phase (resp. quadrature) frequencies. The square (resp. 
round) markers refer to the values computed through numeri- 
cal simulations, of the in-phase (resp. quadrature) solutions. 



the saddle in the full model corresponds to Et — —0.75. 
The Duffing model is an approximation of the motion at 
the bottom of well, i.e. for energies of the order of Eq. For 
higher energies, the Duffing approximation diverges from 
the full model. In particular, the threshold energy in the 
Duffing approximation is Ejj = —1.5, and the in-phase 
solution vanishes for larger energies, as shown in Fig. [2^. 
The reason is that the motion occurs exactly along the 
bissector, where is also located the threshold between the 
wells. When the threshold energy is reached, the particle 
leaves the well. 

To estimate more precisely the domain of validity of 
the Duffing approximation, it is more intuitive to use the 
motion amplitude rather than the energy. Eq. ()13|) gives 
the relation between E and Xq for the l[5c|) and (j8d[) 
solutions: one finds that Xq = 1 corresponds to E = 
— 1.875, which corresponds to the half height of the non- 
approximated well. This means that the interval on which 
the Duffing approximation is realistic is typically —3 < 
E < -1.875. 

To evaluate the stability of this solution, we study the 
motion in its vicinity. We search for solutions of the type: 

X = {Xa + ex) cos {uit + 0^) 
Y = {Xq + sy) cos {uJit + 4)y) 



5 

If this solution is injected into the motion equations, one 
finds that ex+£Y — C where C is a constant. It is easy to 
show that taking C ^ is equivalent to translate from the 
solution (Xq, Xq) to the solution (Xq + ^,Xq + Thus 
we choose ex = —£y = £• The equations of motion, in the 
slow varying envelop approximation (SVEA), become: 

H-^y = [l- 2 ) (14a) 
^ = (14b) 

or: 
with 

The perturbation of the in-phase solution oscillates, and 
thus the in-phase solution is a stable periodic orbit. The 
amplitudes of the motions along X and Y are modu- 
lated in antiphase at the frequency Qi. A phase modu- 
lation, in quadrature with the amplitude modulation, is 
also present, as expressed in (|14b[) . The neighboring so- 
lutions are quasiperiodic solutions with the same main 
frequency oji , and sidebands at the frequency J7i increas- 
ing with the energy (Fig. [^b, solid line). Note that for 
the highest energies {E = —1.5), J7i ~ uji, and thus the 
SVEA is no more valid. However, we already established 
above that this energy domain is not fully relevant, as the 
Duffing approximation is itself too much rough to describe 
the experimental system for such energies. Thus the valid 
domain remains typically —3 < E < —1.875, as discussed 
previously. 

As argued above, the ipx = + solution is the 
same solution as tl)x — i^y when Y is changed in —Y. 
It consists in a particle oscillating on the other bissector. 
This solution {Xq^—Xq) is stable, and neighboring solu- 
tions are an oscillating motion around the main solution. 
The dependence between the frequencies and the energy 
is the same as in the '0x = solution. 

Quadrature solutions 

The last pair of solutions are those where X and Y os- 
cillate in quadrature of phase. The same steps as for the 
V'x = '4'Y lead to the solutions: 

X = Xq cos ((JJ2O (16a) 

Y ^ Xo cos (w2i ± I) (16b) 

with: 

o.|=c.2_(2 + a)^ (17) 

The motion of a particle is a circle around the origin. 
The only difference between both solutions is the rotation 
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direction. The analysis of the motion in the neighboring 
of these solutions leads to an oscillating motion with a 
frequency: 



4w2 V 2 



(18) 



These solutions are stable periodic orbits, and the mo- 
tion in the neighboring is quasiperiodic, with a main fre- 
quency LU2 decreasing with the energy, and sidebands at a 
frequency f22 increasing with the energy (Fig. [21 dashed 
lines). Note that for a circular motion, the kinetic energy 
never vanishes: this is why this motion still exists for total 
energies larger than —1.5, as it does not reach the saddle 
point. 

To summarize the above results, we found that the 
dynamics in the bottom of the main well of the lattice 
is organized around two unstable periodic orbits and four 
stable periodic orbits. As expected, the dynamics in the 
vicinity of these stable periodic orbits follows KAM tori. 
These quasiperiodic regimes can be described as a motion 
at the same frequency as the periodic orbits, perturbed 
by sideband components. Although this system consists 
in two perturbed pendula, the effects of the small per- 
turbation is outstanding: whereas uncoupled pendula can 
oscillate with any relative phase, the coupling introduced 
here allows only four relative phases, even for very small 
motions, where the perturbation tends to zero. To eval- 
uate the pertinence of the above description, let us now 
switch to the results of numerical simulations. 



0.4- 
>< 0.0- 



-0.4- 



>< 



>< 




400 



800 



1200 



0.4 
0.0 
-0.4 




— I — 
10 



20 
T 



30 



40 



0.4 
0.0 



-0.4 



/V /v /' 


^ A " A A \ A 

' i\' A> A 




W' W W > ' 
V' U' U; u; 



(c) 



10 



20 
T 



30 



40 



4 Numerical simulations in the Duffing 
approximation 

The aim of numerical simulations is to determine how the 
previous results evolve far from the bottom of the poten- 
tial wells. Let us first examine the results of numerical 
simulations in the Duffing approximation. The interest of 
such simulations is to give complementary information as 
compared to the previous results, as the shape of the X 
and Y time evolution, or the complexity of the spectra. 

Fig.Oshows the time evolution of X and Y in the vicin- 
ity of the antiphase solution (Fig. [Sb) and of one of the 
quadrature solution (Fig. [3h) for a = 0.5 and E — —2.5. 
In first approximation, it seems that the time evolution in 
both directions follows the same frequency, with a clear 
phase relation between them. However, the small ampli- 
tude modulation (Fig. [3^) shows that other frequencies 
are present. 

More details can be found on the Fourier transform of 
the signals (Fig. [5]). Both dynamics appear to be driven 
by a main frequency shifted as compared to the frequency 
ojQ = y/1 + a at the bottom of the well. The sidebands 
are more than one order of magnitude weaker, and the 
second sidebands are still two orders of magnitude smaller. 
Thus the dynamical regimes along the X and Y directions 
can effectively be approximated as two frequency locked 
oscillations, slightly perturbed by a sideband modulation. 



Figure 3. Motion of the particles as a function of the time 
r — ujot, for a = 0.5 and E — —2.5. (a) large scale motion 
along X; (b) oscillations along X (solid line) and Y (dashed 
line) in the vicinity of the antiphase solution; (c) oscillations 
along X (solid line) and Y (dashed line) in the vicinity of a 
quadrature solution. 



It can be noticed that the spectra of Fig. [5] are asym- 
metric. On Fig. [4^, the high frequency sideband is larger 
than the low frequency one, while it is inverted on Fig.[3)D. 
This sideband imbalance is simply due to the fact that we 
have both amplitude modulation and phase modulation. 
The coupling between these two modulations, expressed in 
(|14a[) and (jl4bp for the in-phasc solution, leads to a high 
frequency component stronger than the low frequency one, 
as observed in Fig. 2^. On the contrary, for the quadrature 
solutions, the high frequency sideband is predicted to be 
weaker than the low frequency one. 

The frequencies obtained by the simulations are re- 
ported on Fig. [21 For low energies, they are identical to 
those predicted theoretically in the previous section. For 
higher energies, small differences appear, in particular for 
the sidebands of the antiphase solution. This originates, 
as discussed previously, from the fact that the SVEA is 
no more valid in this domain, as wi ~ 
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Figure 4. Spectrum of the motion of the particles along the 
X direction, in the vicinity of the (a) anti-phase solution, and 
(b) quadrature solution. These are respectively the FFT of the 
signals of Figs ^jp and (St- Both spectra are represented with 
the same scales, to make the comparison easier. 



5 Behaviors far from the bottom of the wells 

In the above sections, the approximation of the optical lat- 
tice to two coupled Duffing oscillators allows us to show 
that in the vicinity of the bottom of the well, the dynam- 
ics is governed by a synchronization process, very similar 
to the frequency locking. It is well known that this type of 
phenomenon is able to inhibit complex dynamics in dis- 
sipative systems, and if this synchronization applies to a 
large part of the well, it could explain that chaos appears 
only into a very narrow area. We study in this section the 
dynamics of the particles beyond the domain where the 
Duffing approximation is valid. 

Before to examine in details the dynamics of the par- 
ticles in the exact potential, let us remind the global evo- 
lution of the dynamics in the Poincare section. Our phase 

space is 4-dimensional, with directions (^X, Y, X, , but 

because of the energy conservation, the accessible space re- 
duces to a 3D surface. We choose to consider Poincare sec- 
tion at y = with increasing values, and thus, Poincare 

sections are in the 3D space (^X, Y, X^ , and they lie on a 

2D surface Sp, which has the shape looking like a semi- 
ellipsoid. To represent the Poincare sections we could project 
them on the {X, Y) plane, but we use here the more usual 

(^X, X^ plane. It shows the Poincare sections viewed from 

the vertex of the semi-ellipsoid. However , because of the 
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Figure 5. Poincare sections for a = 0.5 and E = —1.02. 



stiff sides of Sp, we have to keep in mind that many curves 
are projected almost at the same location, and thus are 
superimposed. 

For relatively high energies, we know that the Duffing 
approximation is no more valid. Thus it is interesting to 
look if, in spite of that, the behaviors remain those de- 
scribed in the Duffing approximation. As an example, we 
choose arbitrarily to illustrate the following of the paper 
with the behaviors obtained for E — —1.02, i.e. an en- 
ergy well above the domain of validity of the Duffing ap- 
proximation. These behaviors are representative to those 
observed for all large energy values, i.e. typically larger 
than —2. Fig. [5] shows the Poincare section in this situa- 
tion. These results have been obtained through numerical 
resolution of the equations of motion which are derived 
from the potential ([T]), without any approximation. We 
see four distinct domains separated by an X-shaped sepa- 

ratrix. The central point [ X = 0, X = O) corresponds to 



the first trivial solution (|8ap . This solution has been found 
to be unstable, which is compatible with the fact that all 
trajectories move away from it. The second unstable peri- 
odic orbit corresponds to F = 0, and so y = 0, and thus 
it cannot be represented in this Poincare section through 
the Y = plane. 

In the right and left domains, the Poincare sections 
are closed curves around two points with X = 0. As by 
definition of the Poincare section, Y = 0, these points are 
turning points of the motion, and thus correspond to the 
in-phase and antiphase solutions. In the same way, in the 
top and bottom domains, the trajectories are also cycling 
around two points, with X = 0, which correspond to the 
quadrature solutions. Fig. [6] shows the real space trajecto- 
ries of the antiphase solution and of one of the quadrature 
solutions. The antiphase solution appears as a straight line 
on a bissector, along which the particle oscillates with the 
frequency uj'i. The quadrature solution appears as a cir- 
cle around which the particle turns with the frequency 
0^2. The values of iv'i and are reported on Fig. [7^ as a 
function of the energy. As expected, the values predicted 
with the Duffing approximation and those obtained from 
the full model diverge significantly for energies larger than 
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Figure 6. Real space trajectories for the antiphase and the 
quadrature solutions, calculated for a = 0.5, E = —1.02. The 
antiphase solution consists in an oscillation of the particle along 
the straight line, while the quadrature solutions correspond to 
a rotation of the particle along the circle. 



-2. However, the global evolution remains that the main 
frequency decreases as energy is increased. 

Let us now examine the dynamics on the tori close to 
a periodic stable solution. We know that for lower ener- 
gies, it can be considered as periodic oscillations slightly 
perturbed by sidebands. To know if this approximation is 
still valid for larger energies, we look at the spectra of the 
trajectories. Fig. [8] shows the spectrum of two trajecto- 
ries close to the antiphase and quadrature solutions. The 
spectra have definitely the same characteristics as those 
found for the lowest energies: the dynamics is mainly a 
periodic oscillation, with the same frequency in both di- 
rections. This oscillation is perturbed by sidebands, the 
amplitude of which is more than one order of magnitude 
smaller. The second sideband is still one order of magni- 
tude smaller. Thus the behavior in the vicinity of the two 
periodic orbits at large energy remains essentially a peri- 
odic oscillation, where the motion along both directions 
is locked to the same frequency. The frequency beating 
as a function of the energy is reported on Fig. [Tb: as for 
the main frequency, a clear divergence as compared to the 
Dufffing approximation appears at energies larger than -2, 
when the approximation is no more valid. But the global 
evolution remains an increasing of the beating frequency 
with the energy. 

We examine now the behavior of the system in the 
vicinity of the unstable solutions. As discussed in [T7], the 
dynamics in the immediate vicinity of the separatrix is 
chaotic when the energy is large enough. Fig. [TOk shows 
a typical spectrum with the characteristics of chaos, in 
particular a large continuous component. The existence of 
chaotic trajectories in the vicinity of the trivial solutions 
confirms the unstable nature of these solutions. However, 
the chaotic area remains marginal in the lattice we discuss 
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Figure 7. Oscillation frequencies of the particle in the vicinity 
of the periodic orbits. In (a), frequency of the periodic solution; 
in (b), frequency of the beating. The solid (resp. dashed) lines 
are the same as in Fig. O i.e. they refer to the analytical in 
phase (resp. in quadrature) frequencies in the Duffing approx- 
imation. The square (resp. round) markers refer to the values 
computed, through numerical simulations of the full model, for 
the in-phase (resp. in-quadrature) solutions. 



here. For example, in Fig. [51 the chaotic area is so small 
that it cannot be visible at the scale of the representation. 

Fig. [9] shows the temporal dynamics on tori passing 
close to the {X = 0, X = 0) point, just outside the chaotic 
area, on the antiphase solution side in (a), and on the 
quadrature side in (b). As expected in the vicinity of an 
unstable solution, these dynamics do not seem to be linked 
to the solution itself. On the contrary, they appear as an 
evolution of the regimes generated by the corresponding 
stable periodic orbit. Indeed, they have the main char- 
acteristics of the regimes observed in the vicinity of the 
periodic cycles. In particular, the two regimes differ by 
the phase difference between the oscillations in X and Y, 
and the amplitude modulations for the motion in X and 
in Y are always in antiphase, i.e. the particle oscillates for 
a while along the X-axis and then moves to the F-axis, 
and vice- versa. Fig. [TOb andfTOb shows the spectra associ- 
ated with these dynamics. The dynamics in the X and Y 
directions are still synchronized, with one main frequency 
and sidebands well below this main frequency. However, 
some differences appear as compared to the vicinity of the 
stable solutions: in particular, the amplitude of the side- 
bands is larger, with less than one order of magnitude as 



D. Hennequin and P. Verkerk: Synchronization in non dissipative optical lattices 



9 




10 1 • 1 • 1 • 1 

0.6 0.8 1.0 1.2 



Figure 8. Spectrum of the motion along the X direction, in 
the vicinity of the (a) antiphase solution, and (b) quadrature 
solution for a = 0.5, E — —1.02. The corresponding trajecto- 
ries give the Poincare sections which are the first circles around 
the periodic solutions in Fig.O In (a), Poincare section crossing 
the X = axis in X = 1.46. In (b), Poincare section crossing 
the X = axis in X = 1.195. 



compared to the main frequency, and several harmonics 
have a non negligible amplitude. 

We have now a good picture of the motion of an atom 
in a well of our lattice. Whatever its energy is, its mo- 
tion is mainly a periodic oscillation distorted by a slow 
drift. This drift vanishes on the periodic stable solutions, 
and increases when one moves away from these solutions. 
The frequency of the oscillation, together with that of the 
sidebands, evolve with the energy, but evolve also in the 
phase space of a given energy. Fig. [TT] shows the evolution 
of the behavior when one moves away from the stable solu- 
tions, for a = 0.5 and E = —1.02. On the whole accessible 
phase space, the dynamics is still essentially an oscillation 
with the same frequency in both directions, perturbed by 
sidebands. Fig. [TT] shows that the main frequency evolves 
continuously between the resonance frequencies of the dif- 
ferent stable solutions. In particular, for the four domains 
delimited by the separatrix, the main frequency tends to 
that of the unstable solution for trajectories at the edge of 
the domain. The beating frequency also evolves monoton- 
ically in each domain, tending to zero when the separatrix 
is approached. So the transition between each periodic so- 
lution occurs without any discontinuity. 




T ' 1 ^ 1 ' 1 



100 200 300 

T 

"lIMilMI 

T ' 1 ' 1 ' 1 

100 200 300 

T 

Figure 9. Motion of the particles as a function of the time 
r = Luot, for a — 0.5 and E = —1.02, for initial conditions 
close to the unstable periodic orbits. In both plots, X and 

Y trajectories are centered on 0, but for sake of clarity, the 

Y trajectories have been shifted by vr, so that they appear 
to be centered on tt. The corresponding trajectories are on 
tori centered on (a) the antiphase periodic orbit and (b) a 
quadrature periodic orbit. 

6 Conclusion 

We have shown in [17] that the nature and the complex- 
ity of the motion of an atom in an optical lattice - or 
similarly of a classical particle in the corresponding exter- 
nal potential - can change drastically as a function of the 
lattice properties. We study the fundamental mechanisms 
leading to so important differences in the dynamics, and 
especially the absence of chaos in the potential wells of a 
red detuned square optical lattice. We adopt an approach 
quite unusual in the domain of conservative chaos, because 
our aim is to find experimental tools able to characterize 
more precisely the complex conservative dynamics, and 
in particular able to distinguish between different com- 
plex behaviors. In square lattices, atoms traveling between 
sites follow two different chaotic behavior, depending on 
the lattice blue or red detuning. This difference proba- 
bly originates in the behavior of atoms inside the wells: 
full chaos appears when the laser frequencies are blue de- 
tuned, whereas chaotic trajectories are quasi inexistent 
when they are red detuned. The former is not surprising: 
the edges of a well are far from harmonicity, and the cou- 
pling between the two directions, together with the non- 
linearities, become very high. On the contrary, the mecha- 
nism leading to the absence of chaos in the latter needed to 
be clarified. We show here that synchronization between 
the motion in the two directions of space inhibits chaos, 
through a mechanism very similar to that of frequency 
locking in dissipative systems. Indeed, the frequencies in 
both directions of space are degenerate at the bottom of 



10 



D. Hennequin and P. Verkerk: Synchronization in non dissipative optical lattices 




0.5 0.6 0.7 0.8 0.9 1.0 
CO/COq 

Figure 10. Spectrum of particle trajectories close to the un- 
stable periodic cycles, in (a) the chaotic area, (b) on the side 
of the antiphase solution and (c) on the side of a quadrature 
solution. 
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Figure 11. Evolution of the main frequency cj (squares) and 
the beating O (crosses) for different trajectories, for a — 0.5 
and E = —1.02. Trajectories are identified through the point 
where they intersect a line linking the quadrature and the an- 
tiphase solution. D is the distance between the intersection 
point and the separatrix. The quadrature (resp. antiphase) so- 
lution is in D = -0.4 (resp. D = 0.6). 



the well. For atoms with higher energies, the degeneracy 
should be broken, but the coupling overcomes the anhar- 
monicity and leads to a motion with one single frequency. 
Synchronization is not strictly frequency locking, as our 
system is conservative, and there is no dissipation to com- 
pensate for frequency pulhng. So the periodic oscillations 
are modulated in phase and in amplitude, i.e. the main 
frequency comes along with small sidebands. In the phase 
space, the amplitude modulation can increase, and even 
reach 100% in the vicinity of the unstable periodic orbits, 
but the description in terms of frequency locking remains 
valid on most of the phase space. Finally, in our conserva- 
tive lattice, synchronization appears to be at the origin of 
the absence of chaos, following the same mechanisms as in 
dissipative systems. This intrinsic origin of this inhibition 
is not a local behavior, i.e. the lack of chaos will remain 
whatever the lattice parameters are, except for the lattice 
frequency. In particular, the quantum dynamics should be 
studied in another lattice, as e.g. the blue detuned square 
lattice, where chaos is fully developed. 

This analysis shows that the dynamics of conserva- 
tive systems can be characterized using deterministic tools 
rather than the traditional statistical approaches. Atoms 
in optical lattices are a good model system to apply these 
tools and test them on experiments. Future studies could 
refine the present results, in particular by checking that 
this approach is still valid when, at the main resonance, 
the frequencies in both spatial directions are not equal, 
but only multiple. This can be achieved with the blue de- 
tuned square lattice, by choosing adequately a. 
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